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The behaviour of the topological susceptibility in QCD with two colours and 8 flavours of quarks is studied 
at nonzero temperature on the lattice across the finite density transition. It is shown that its signal drops at a 
(pseudo-) critical chemical potential fi c . The Polyakov loop and the chiral condensate undergo their transitions 
at the same value. Pauli blocking supervenes at a value of the chemical potential larger than \x c . 



1. INTRODUCTION 

The phase space of QCD in the temperature 
and quark chemical potential fx plane is shown 
schematically in Fig. ^ Many of the transi- 
tion lines and phases displayed in the Figure are 
theoretical predictions with little experimental 
verification. We have studied the transition as 
the chemical potential is varied (dashed arrow in 
Fig.HJ at a fixed temperature T that is most likely 
above the region where a diquark condenses. 

We wanted to understand the fate of the topo- 
logical susceptibility x across the transition and 
to decide whether its possible change occurs at 
the same value of y, where the vector chiral sym- 
metry is restored and the Polyakov loop signal 
rises £Q. 

It is known that at the finite temperature tran- 
sition the signal of x drops abruptly for Yang- 
Mills theory 2 3 as well as for QCD with sev- 
eral values of the flavour number Nf |3|. The 
present is the first work where a similar study is 
performed at the finite density transition (the full 
account can be found in Ref. 

We have performed the study on the lattice 
by simulating the theory with two colours in or- 
der to avoid the sign problem that in the three 
colour theory makes the importance sampling 
method impracticable (for methods to overcome 
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Figure 1. Sketch of the T-fi phase diagram in 
QCD with two colours. We have studied the tran- 
sition crossed by the dashed arrow at finite tem- 
perature as the density is varied. 



this problem see [H] and references therein). We 
expect that this modification of true QCD has lit- 
tle impact on the results concerning gluon prop- 
erties and in particular instanton physics. The 
main differences between two and three colour 
QCD regard aspects of the theory that we have 
not studied, like the hadron spectrum and the 
diquark condensed phase |7l8j . In fact in the 
two colour theory mesons are degenerate with 
baryons (which are made of two quarks) and the 
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diquark condensate phase becomes superfluid in- 
stead of superconducting because two quarks in 
two colour QCD can create a colour blind con- 
densate. 

2. SIMULATION AND OBSERVABLES 

We have simulated the model on a 14 3 x 6 
lattice at an inverse bare lattice gauge coupling 
(3 = A/g 2 = 1.5 and quark mass am = 0.07 in 
units of lattice spacing a. The Wilson action was 
used for gauge fields and the standard action for 
8 flavours of quarks in the staggered formulation. 

The Hybrid Molecular Dynamics algorithm 
was chosen to update configurations. They we- 
re separated typically by 50-100 steps of the al- 
gorithm in order to well decorrelate the topol- 
ogy |9I10I11| . We followed the behaviour of sev- 
eral observables across the transition shown in 
Fig. CJ chiral condensate (tpip), Polyakov loop or 
Wilson line P, topological susceptibility x> bary- 
onic density pb and average plaquette Tr d/2 

2.1. Topology 

The measurement of the topological suscepti- 
bility requires a careful treatment of composite 
operators. This quantity is defined as 

X = J d A x{T{Q(x)Q(0)}) , (1) 

where Q(x) is the topological charge density oper- 
ator that appears in the r.h.s. of the axial singlet 
anomaly, x 1S related to the 77' mass through the 
Witten-Veneziano mechanism |12I1H| and to the 
chiral condensate in the massless limit |14| . 

Eq.Q can be rewritten in a more compact 
form, 



where the total topological charge Q is defined 
as Q = J d 4 xQ(x) and V is the spacetime vol- 
ume. This expression has only a formal meaning 
because it must be supplemented with a multi- 
plicative renormalization for each power of the 
topological charge operator Q |15I16| as well as a 
careful treatment of the contact divergences that 
appear when the two operators are evaluated at 
the same spacetime point |17j . 



A detailed description of the procedure that we 
have followed to extract \ can be found in 0. 
In a nutshell, one has to define a lattice regular- 
ization Ql{x) of the topological charge density 
operator. The corresponding total lattice topo- 
logical charge is Ql = Ylx Ql(x) and the lattice 
topological susceptibility is xl = (Ql)/V- This 
quantity is related to the physical susceptibility 
X by the general expression |15I16I17| 

XL = Z 2 a 4 X + M , (3) 

where Z is a multiplicative renormalization of Ql , 
a is the lattice spacing and M is an additive 
renormalization constant. In presence of com- 
posite operators the usual renormalization pro- 
gram in a renormalizablc theory is not enough to 
fully renormalize a Green's function. This is the 
origin of Z (even in the case of the pure gauge 
theory where the topological charge operator has 
no anomalous dimensions). On the other hand 
part of the contact divergences that appear in the 
product of the two topological charge operators 
must be subtracted. This subtraction is defined 
as M = xl\q=o, i- e - M is the value of xl in 
the sector of zero topological charge. This con- 
dition guarantees the obvious requirement that x 
vanishes in that sector. 

2.2. Calculation of Z and M 

One has to calculate the two constants Z and 
M and insert them into Eq.® in order to ex- 
tract x- We shall present the results in the form 
of the ratio x(/ i )/x(/ i = 0) as a function of \x. 
Therefore we need not calculate Z and only M 
is required. This additive renormalization con- 
stant was computed by using a non-perturbative 
method |18ll9lri| . It is enough to calculate it for 
one single value of the chemical potential because 
M is independent of infrared properties. 

3. RESULTS 

In Fig. |3 we show the resulting behaviour of 
the topological susceptibility as a function of the 
chemical potential. A clear cut drop is seen at the 
position a/i c = 0.175(5). This sudden fall is an 
indication of the effective restoration of the axial 
singlet symmetry |2U| . 
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Figure 2. The ratio x(/i)/x(/i = 0) as a function 
of a/i showing the abrupt drop of the signal of the 
topological susceptibility at a (pseudo-) critical 
value of the chemical potential. 



We do not show the analogous figures for the 
Polyakov loop and chiral condensate behaviours 
(see Ref. jS])- Instead in Fig. the derivatives 
of the three quantities (topological susceptibil- 
ity, Polyakov loop and chiral condensate) are su- 
perimposed to give evidence of the coincidence 
of the three transitions: they all happen at the 
same (pseudo-)critical value a/i c = 0.175(5). Sin- 
gle points are the result of the discrete deriva- 
tives computed from the measured data while 
the lines are the derivatives of the interpolations 
made with a natural cubic spline on each data 
set. Since the theory deconfines at the same \i 
where (ipt/)) vanishes, the path in the T-fj- dia- 
gram followed in our study possibly lies above 
the superfiuid phase (as indicated in Fig.^l |21|. 

3.1. Numerical results 

In order to give an estimate of the main re- 
sults in physical units, we have done a second 
simulation of the model at zero chemical poten- 
tial by varying the inverse bare lattice coupling 
P on the same lattice volume and have measured 
the Polyakov loop and chiral condensate. The 
data clearly reflected the existence of a critical 
j3 c — 1.594(6). It was then assumed that the 
corresponding transition temperature T c lies in 
the range 100-200 MeV. Hence a value for the 
lattice spacing at j3 c was obtained. This value 
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Figure 3. Derivatives with respect to the chemi- 
cal potential fx of the normalized topological sus- 
ceptibility (circles, continuous line and left ver- 
tical axis), Polyakov loop (up triangles, dotted 
line and right vertical axis) and chiral conden- 
sate (down triangles, dashed line and right ver- 
tical axis). The two vertical axes are expressed 
in units of the lattice spacing a. The three sets 
of data have been slightly shifted horizontally to 
avoid the overlapping of various symbols and er- 
ror bars. Lines are the result of a spline interpo- 
lation. 



was then run, by making use of the two-loop 
beta function, from f3 c to our working [3 — 1.5 
thus obtaining the lattice spacing at this beta 
a(/3 = 1.5) = 0.64(4)(±ig) fm (errors derived from 
the imprecision on j3 c and from the inaccuracy of 
the estimate for T c respectively). 

This result allowed us to assign physical units 
to the various parameters of our simulation: fie = 
54(2) (4) (18) MeV (errors are respectively due 
to the estimate of a^ c , the imprecision on (3 C 
and the rough approximation of T c ) and T = 
51(4) (17) MeV (error on f3 c and on T c respec- 
tively). 

Notice that the numerical value of /z c is roughly 
compatible with half of the mass of the lightest 
baryon (in the two colour theory it is made of 
two quarks and degenerated with the pion) . How- 
ever we take the above numbers with caution be- 
cause of the large systematic errors derived from 
the small and coarse lattice volume, the (wrong) 
gauge group and the inexact updating algorithm. 
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3.2. Pauli blocking 

By repeating the simulation at very high val- 
ues of the chemical potential, in principle one 
could test the perturbative region of the phase 
diagram. However the advent of the so-called 
Pauli blocking makes that not viable. As fi in- 
creases, more and more fermions are placed in 
the lattice volume. Since it contains a count- 
able number of sites, the Pauli principle imposes 
a maximum allowed number of fermions that can 
be placed on them. In our case the maximum 
allowed density is one baryon per site [22] • This 
maximum is attained for et/i s ks 1.2. After that 
point, fermions are frozen and the theory becomes 
entirely quenched. We have verified this state- 
ment by calculating the baryon density and the 
average plaquette. Their data are shown in Fig.Q] 
Pb becomes 1 and stays constant for all /i > fj, s . 
Also the average plaquette stays constant after 
/x s and this constant coincides with the value ob- 
tained from a separate Monte Carlo study where 
the pure gauge theory was simulated (represented 
by the down triangle in Fig. E}. We stress that 
Pauli blocking is a lattice effect with no counter- 
part in the continuum. 



REFERENCES 



1 



0.75 




0.25 



Figure 4. Baryon density per space site and av- 
erage plaquette as a function of ap. The down 
triangle is the result of a separate simulation for 
the pure gauge theory. 
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